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Abstract 

We present the exact analytical solution of the radial Schrodinger equation for the deformed 
Hulthen and the Morse potentials within the framework of the Asymptotic Iteration Method. The 
bound state energy eigenvalues and corresponding wave functions are obtained explicitly. Our 
results are in excellent agreement with the findings of the other methods. 
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I. INTRODUCTION 



In recent years, the energy eigenvalues and corresponding eigenfunctions between interac- 
tion systems have raised a great deal of interest in relativistic quantum mechanics as well as 
in non-relativistic quantum mechanics. The exact solution of the wave equations (relativis- 
tic or non-relativistic) is very important since the wave function contains all the necessary 
information regarding the quantum system under consideration. Analytical methods such as 
the super-symmetry (SUSY) [l] and the Nikiforov-Uvarov method (NU) [2| have been used 
to solve wave equations exactly or quasi-exactly for a given potential and these potentials 
are in general either the polynomial type like the Coulomb, the harmonic oscillator 3| or 
the exponential type such as the Eckart or Hulthen j^, 4, (|, Morse or a combination of 
these potentials. 

The Hulthen and Morse potentials we consider in this paper have a wide range appli- 
cations. The Hulthen potential has a very important role in describing the bound states 
or continuum states between interaction particles in the relativistic quantum mechanics as 



well as the non-relativistic quantum mechanics. The Hulthen potential has been solved for 
;he bound states with the NU 

a , the SUSY 0, y, y, Q 

and the shifted 1/N expansion 
11 | methods. It is also solved for the relativistic wave equations such as Dirac 12. T13I. 
Klein-Gordon (KG) 3, Q and the Duffin-Kemmer Petiau (DKP) [17] equations. 

The Morse potential has raised a great deal of interest over the years and has been one 
of the most useful models to describe the interaction between two atoms in a diatomic 
molecule. In this respect, Morse potential has been solved by the super-symmetry (SUSY) 
[l|, the Nikiforov-Uvarov method (NU)[2], the hypervirial perturbation method (HV) 18], 
the shifted and modified shifted 1/N expansion methods 19| as well as the variational 



method 



20|. 



Since these potentials have been extensively used to describe the bound and the con- 
tinuum states of the interactions systems, it would be interesting and important to solve 
the non-relativistic radial Schrodinger equation for the Hulthen and Morse potentials. Re- 
cently, an alternative method called as the Asymptotic Iteration Method (AIM) for solving 
second-order homogeneous linear differential equation has been developed by Ciftgi et al. 



21 



22j and has been applied to solve the non-relativistic radial Schrodinger equation or the 



relativistic Dirac equation. 
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In this paper, our aim is to solve the deformed Hulthen potential and the Morse potentials 
to obtain the energy eigenvalues and corresponding eigenfunctions within framework of the 
AIM. In the next section, the asymptotic iteration method (AIM) is introduced. Then, 
in section IIII[ the Schrodinger equation is solved by the asymptotic iteration method for 
the deformed Hulthen potential and the Morse potential. The AIM results are compared 
with the findings of the other methods. Finally, section IIVI is devoted to the summary and 
conclusion. 



II. OVERVIEW THE ASYMPTOTIC ITERATION METHOD (AIM) 



A. Energy Eigenvalues 



AIM is briefly outlined here and the details can be found in references 2jJ, |22j, |23j. AIM 
is proposed to solve the second-order differential equations of the form 



y" = X (x)y' + s (x)y 



(1) 



where Xq(x) ^ and the prime denotes the derivative with respect to x. The variables, 
so(x) and Xq(x), are sufficiently different iable. The differential equation (CQ) has a general 
solution [2l|] 



y(x) = exp 



/x \ r px / pxx \ 

a(xi)dxi J C*2 + C\ l exp I / \Xq{x2) + 2a{x2)}dx2 ) dx\ 



(2) 



for sufficiently large k, k > 0, if 

s k (x) s k -i(x) 



X k (x) Afe_i(x) 



a(x), k = 1, 2, 3, 



(3) 



where 



X k (x) = X'^x) + s k ^{x) + A (a;)A fc _i(a;) 

s k {x) = + s o(x)X k -i{x), k = 1,2,3,. 



(4) 



Note that one can also start the recurrence relations from k = with the initial conditions 



A_i = 1 and s_i = 24j. For a given potential, the radial Schrodinger equation is converted 
to the form of equation ([1]). Then, Sq(x) and Xq(x) are determined and s k (x) and X k (x) 
parameters are calculated by the recurrence relations given by equation (jl]). 



The termination condition of the method in equation ([3]) can be arranged as 



A fe (x) = A fe (x)s fe _i(x) - \ k -i(x)s k (x) = k = 1,2,3,... (5) 

The energy eigenvalues are obtained from the roots of the equation §5$) if the problem 
is exactly solvable. If not, for a specific n principal quantum number, we choose a suitable 
xq point, determined generally as the maximum value of the asymptotic wave function or 
the minimum value of the potential 2l|, |24j, and the approximate energy eigenvalues are 



obtained from the roots of this equation for sufficiently great values of k with iteration. 



B. Energy Eigenfunctions 

In this study, we seek the exact solution of radial Schrodinger equation for which the 
relevant second order homogenous linear differential equation takes the following general 
form [22|, 

s „ / ax N+1 (m + l)\ ,. . wx N , s . . 



the following general formula for the exact solutions for y n (x) is given by [22J, 

y n {x) = (-1)" C 2 (N + 2) n (a) n 2 F l (-n, p + n;a; bx N+2 ) (7) 

whprp (rr) -JM rr - j _ (2m+l)b+2a 

wnere [a) n - r((r) , a - N+2 ana p - (A r +2 )6 ■ 



III. CALCULATION OF THE ENERGY EIGENVALUES AND EIGENFUNC- 
TIONS 

The motion of a particle with the mass m is described by the following Schrodinger 
equation: 



-h 2 ( d 2 



2 d 



2m V dr 2 r dr r 2 



1 g 
sin 6 1 9^ 



1 

sin 2 (90 



+ ^(r)U n£m (r,0, 



(8) 



The terms in the square brackets with the overall minus sign is the dimensionless angular 



momentum squared operator, L 2 . Defining ^ n e m (r, 9, 



u n t(r)Y em ( 



we obtain the 
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radial part of the Schrodinger equation: 



d 2 2d 
dr 2 r dr 



2m 



V(r) + 



2mr 2 
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2mE 

Untir) + —r^-u ni {r) = 



It is sometimes convenient to define u n i(r) and the effective potential as follows 



Rni(r) 



eff 



V(r) + 



£+i)h 2 

2mr 2 



(9) 



(10) 



Since 



d^ + 2d_\ R nl {r) 
dr 2 r dr I r 



r dr 2 



RnAr) 



The radial Schrodinger equation given by equation (JOJ) follows that 

d 2 R n£ (r) 2m r 



(11) 



(12) 



Instead of solving the partial differential equation ([8]) in three variables r, 9 and <p, we 
now solve a differential equation involving only the variable r and the angular momentum 
parameter i = (s-state solution). 



A. Deformed Hulthen Potential Case 

The deformed Hulthen potential j^] is defined as 



e 



-Sr 



V DH (r) = -Ze 2 5 T -^ (13) 

In this formula, Z, 5 and q are respectively the atomic number, the screening parameter and 
the deformation parameter determining the range for the deformed Hulthen potential. The 
deformed Hulthen potential reduces to the Hulthen potential form for q — 1, to the standard 
Wood-Saxon potential for q — — 1 and to the exponential potential for q = 0. 

Inserting the Hulthen potential given by equation ( TlBl) into equation [12] and using follow- 
ing ansatzs 

_ 2mE 2m^Z 

h 2 s 2 ' p h 2 8 ' y ' 

and we apply a transformation to 5r = x to the deformed Hulthen potential. The radial 
Schrodinger equation takes the following form: 

d 2 R(x) 



^ v- e! + 3 V^) w = (15) 



5 



If we rewrite equation [15] by using a new variable of the form y = e x , we obtain 



d 2 R{y) , ldR(y) 



dy 2 



+ - 



V dy 



+ 



R{y) = 



(16) 



v* y(i-sy). 

In order to solve this equation with AIM, we should transform this equation to the form of 
equation [TJ Therefore, the reasonable physical wave function we propose is as follows 



R(y) = y%l - qy)f(y) 



(17) 



If we insert this wave function into the equation [16], we have the second-order homogeneous 
linear differential equations in the following form 



d 2 f(y) 
dy 2 



(2eq + 3q)y-(2s + l) 



dm 

dy 



+ 



2eq + q - (3 2 

- qy) 



f(y) 



(18) 



1/(1 - qy) 

which is now amenable to an AIM solution. By comparing this equation with equation [1] 
we can write the Xo(y) and so(y) values and by means of equation 0], we may calculate \k(y) 
and Sk(y). This gives: 

,'2eqy + 3qy-2e-l 

s 
Ax 



2/(1 - qy) 

2eg + g-(3 2 

y(i - ??/) 

2 + /? 2 gy 2 - 8 e 2 ?/g + 4 £ Vg 2 - y/5 2 + 12gVe + 11 gV -18eyg + 6e-7g?/ + 4£ 



y 2 {-l + qyY 



(2eq + q-p 2 ) (-2 + 5qy-2e + 2eyq) 



qy)' 



(19) 



Combining these results with the quantization condition given by equation [5] yields 

lq-(3 2 



So 

Ao 
£i 
Ai 

A 2 

etc 



£i 
Ai 

A 2 

A 3 



£2 



2 g 




14g- 




4 g 




19g- 




6 g 



(20) 



When the above expressions are generalized, the eigenvalues turn out as 

1 / q(n + l) 2 - (3 2 ' 



(2n + 2) 



q 



n = 0,1,2,3, 



(21) 
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Using equation HH we obtain the energy eigenvalues E n , 

2 



2m 



me-'Z 



(n + l)8 



n = 0,1,2,3, 



[n + l)q 2 

In the atomic units (h — m — e—1) and for Z = 1, equation 1221 turns out to be 



(22) 



E n = — - 



1 

nq 



n5 



n = n + 1 n = 1, 2, 3, ... 



(23) 



It may thus be seen that the energy eigenvalue equation is easily obtained by using AIM. 
In order to test the accuracy of equation [231 we calculate the energy eigenvalues for Z = 1, 
(7 = 1 and several values of the screening parameter. The AIM results are compared with 

n n 

the Nikiforov-Uvarov |6j and the shifted 1/N [ll|] expansion methods in Table [B It may 
here be seen that the AIM results are in excellent agrement with the findings of the other 
methods. 

Now, as indicated in Section [Til we can determine the corresponding wave functions by 
using equation [3 When we compare equation [6] and equation [T8| we find iV = —1, b = q, 
a = q, and m = ^f^- Therefore, we obtain p = 2e + 2 and a = 2e + 1. Having determined 
these parameters, we can easily find the eigenf unctions f n (y) by using equation [7| as follows 

f n (y) = (-ir r( ^" + n + 1) 2 F 1 (-n, 2e n + n + 2; 2e n + 1; qy) (24) 



Finally, we can write the total radial wave function as below, 



where C„ 



-1 



R n = NC n y £n (l - qy) 2 F 1 {-n, 2e n + n + 2; 2e n + 1; qy) 



r(2e„+n+l) 
1 r(2 £n +l) ' 



(25) 



B. Morse Potential Case 



The Morse potential is defined as 



V M orse{r) = D e (e- 2ax - 2e~ ax ) 



(26) 



with x = (r — r e )/r e and a = ar e . Here, D e and a denote the dissociation energy and Morse 
parameter, respectively. r e is the equilibrium distance (bound length) between nuclei and a 
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is a parameter to control the width of the potential well. Inserting equation [26] into equation 
[121 and using the following ansatzs 

2 _2jMiE _ 2yr 2 e D e 

£ ~ h 2 ' P ~ h 2 1 ] 

The radial Schrodinger equation takes the following form: 

cPR„(x 



+ (-e 2 + 23 2 e~ ax - 3 2 e~ 2ax ) R n (x) = (28) 



dx 2 

If we rewrite equation [28] by using a new variable of the form y = e~ ax , we obtain 



d 2 R n {y) , ldR n (y) 



e 2 1 2f3 2 1 _ 3 2 
a 2 y 2 a 2 y a 2 



R n {y) = o (29) 



dy 2 y dy 

In order to solve this equation with AIM for I = 0, we should transform this equation to 
the form of equation [T] Therefore, the reasonable physical wave function we propose is as 
follows 

R n {y) = yte-^My) (30) 

If we insert this wave function into the equation [29], we have the second-order homogeneous 
linear differential equations in the following form 

d 2 f n (y) = r Wy-2e-a \ df n (y) + ( 2e(3 + a(3 - 2(3 2 \ 
dy 2 V a y J dy \ ya 2 J 

which is now amenable to an AIM solution. By comparing this equation with equation [H 
we can write the Xo(y) and s (y) values and by means of equation |H we may calculate Xk{y) 
and Sk(y). This gives (the subscripts are omitted): 

V a y 

f2ep + ap-2p 2 
s ° = 2 

V ya 2 

-3 (3 ay + 6e a + 2 a 2 - Qf3ye - 2 f3 2 y + Ae 2 + Af3 2 y 2 

a z y z 

8 (2e + a-2{3)(-a-E + l3y) 

si — i ^ [oz 

a A y l 

. . . etc 
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Combining these results with the quantization condition given by equation [5] yields 



so 

Ao " 


Sl 

Ai 


» e = 




Sl 

Ai = 


S2 

A 2 


» ei = 


3a 

- t +/j 


A 2 ~ 


S3 

A 3 


» £2 = 


5a _ 

—+f> 


etc 









(33) 



When the above expressions are generalized, the eigenvalues turn out as 

e n = /5-(n + i)a n = 0, 1, 2, 3, ... (34) 



Using equation [271 we obtain the energy eigenvalues E 



^ 2 r 



n 



(3- (n+ ~)a 



(35) 



2/ir e 2 

we calculate the energy eigenvalues of the H 2 diatomic molecule. The AIM results are 
compared with those obtained by the hypervirial perturbation method (HV) [18], the shifted 
1/N and modified shifted 1/N expansion methods [19J for the H 2 diatomic molecule in Table 
im As it can be seen from the results presented in these tables that the AIM results are in 
good agreement with the findings of the other methods. 

Now, As indicated in Section [Til we can determine corresponding wave function by using 
equation [7J When we compare equation [6] and equation [3T| we find N = — 1, a = -, and 
m = %p2. Therefore, we find a = — + 1. For b — > we can take the limit in equation [7J 
using the limit relation 

Hm 2 F 1 (—n, \ + a; c; zb) = iF x (— n; c; z) (36) 
fe^o b 



Consequently, the solution of equation [31] can easily find 



f(\ < ^n T C J t + n + 1 ) 2e„ 2/3 



Thus, we can write the total radial wave function as below, 



R n {y) = Ny^e-^F.i-n, ^ + 1; (38) 

a a 



where iV is normalization constant. 



9 



IV. CONCLUSION 



We have shown an alternative method to obtain the energy eigenvalues and corresponding 
eigenfunctions of the deformed Hulthen and the Morse potentials within the framework of 
the asymptotic iteration method. We have calculated the energy eigenvalues for Hulthen 
potential with Z = 1, q = 1 and several values of the screening parameter. The AIM results 
are compared with the Nikiforov-Uvarov j(| and the shifted 1/N expansion methods in 
Table [H Furthermore, We have calculated the energy eigenvalues for H2 diatomic molecule 



and compared the hypervirial perturbation method 18( , the shifted 1 /N and modified shifted 
1/N expansion methods [3] in Table [III As it can be seen from the results presented in 
these tables, the AIM results are in good agreement with the findings of the other methods. 

The advantage of the asymptotic iteration method is that it gives the eigenvalues directly 
by transforming the radial Schrodinger equation into a form of y" =\o(r)y' + So(r)y. The 
method presented in this study is a systematic one and it is very efficient and practical. It 
is worth extending this method to the solution of other interaction problems. 
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n -E n (l/N) [11] -E n (NU)\Q\ —E n (AIM) 



5 = 0.002 



1 


0.4990005 


0.4990005 


0.4990005 


2 


0.1240020 


0.1240020 


0.1240020 


3 


0.0545601 


0.0545601 


0.0545601 


4 


0.0302580 


0.0302580 


0.0302580 


5 




0.0012500 


0.0012500 






5 = 0.01 




1 


0.4950125 


0.4950125 


0.4950125 


2 


0.1200500 


0.1200500 


0.1200500 


3 


0.0506681 


0.0506681 


0.0506681 


4 


0.0264501 


0.0264500 


0.0264500 


5 


0.0153128 


0.0153125 


0.0153125 






5 = 0.05 




1 


0.4753125 


0.4753125 


0.4753125 


2 


0.1012503 


0.1012500 


0.1012500 


3 


0.0333746 


0.0333681 


0.0333681 


4 


0.0113035 


0.0112500 


0.0112500 


5 




0.0028125 


0.0028125 



1 

2 
3 
4 
5 



0.4049962 
0.0450856 



6 = 0.2 
0.4050000 
0.0450000 
0.0005556 
0.0112500 
0.0450000 



0.4050000 
0.0450000 
0.0005556 
0.0112500 
0.0450000 



n 

TABLE I: The comparison of the AIM results (present work) with the findings of the 1/N [ll| and 
NU [6i] methods for the s-states energy eigenvalues of the Hulthen potential for several values of 
screening parameter 5. ^2 



n £ AIM HV [18] Modified Shifted 1/N [19] Shifted 1/N [19] 
-4.47601 -4.47601 -4.4760 -4.4749 

5 -2.22052 -2.22051 -2.2205 -2.2038 

7 -1.53744 -1.53743 -1.5374 -1.5168 



TABLE II: For the H2 diatomic molecule, the comparison of the energy eigenvalues (in eV) 
obtained by using AIM with other methods for different values of n. Potential parameters are 
D = 4.7U6eV, a = 1.9425(A )- 1 , r e = 0.7416A , he = 1973.29eI/^° and pt = 0.50391amu. 
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